Counting Statistics of Many-Particle Quantum Walks 
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We study quantum walks of many non-interacting particles on a beam splitter array as a paradig- 
matic testing ground for the competition of single- and many-particle interference in a multi-mode 
system. We derive a general expression for multi-mode particle-number correlation functions, valid 
for bosons and fermions, and infer pronounced signatures of many-particle interferences in the count- 
ing statistics. 
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I. INTRODUCTION 



Quantum walks (QWs) provide the natural - yet uni- 
tary - extension of the classical random walk (RW) to the 
quantum regime JTJ [2J, where the superposition prin- 
ciple allows a wave function to explore different paths 
simultaneously. The resulting interference in the transi- 
tion amplitudes is cause to effects like the linear increase 
in time of the mean travelled distance of the quantum 
walker 0|4). This ballistic behavior contrasts with the 
diffusive behavior of a classical particle performing a 
RW, whose mean travelled distance increases with the 
square root of time. Single-particle QWs, as realized in 
various experiments (5H8||, can also be mimicked with 
classical wave mechanics 119 H1211 . In contrast, QWs with 
two particles [13-18] show purely quantum features due 
to many-particle interference or entanglement that no 
classical system can provide. For larger particle num- 
bers, it is unknown how the emerging hierarchy of many- 
particle interference effects in terms of the number of 
contributing particles affects physical observables. We 
will derive a general expression for multi-mode parti- 
cle number correlation functions which gives access to 
this hierarchy, and provide first examples as manifest 
in prominent interference structures in two-mode num- 
ber correlation functions. While the QW constitutes a 
paradigmatic testing ground for the study of such many- 
particle interference, the formalism we present applies 
to any linear, non-interacting scattering scenario. 

In a many-particle QW, the distinction between ef- 
fects of single-particle interference and those due to gen- 
uine many-particle interference is crucial. Since distin- 
guishable particles show no many-particle interference, 
a scenario with non-interacting, distinguishable parti- 
cles Jl9 1 is equivalent to many subsequent, independent 
single-particle walks. Given the position probability dis- 
tribution of an individual particle, standard combina- 
torics suffices to derive the many-body counting statis- 
tics. Therefore, the differences in the behavior of distin- 
guishable and indistinguishable particles highlight true 



multi-particle interference. We will see that these differ- 
ences are not detectable by the mean particle number, 
but only by observables which capture many-particle 
correlations. The counting statistics thus encodes an 
unambiguous signature of many-particle coherence and 
interference and is readily accessible in the experiment. 



II. FRAMEWORK 

We consider a many-particle QW [20] in one dimen- 
sion, physically realized by N particles which propagate 
through an array of L unbiased beam splitters (BS) as 
depicted in Fig.[T][21 , 22]. The transverse position of the 
beam splitters in the array defines the particle position, 
and subsequent BS rows represent the discrete steps in 
time. A particle that passes a single BS will leave it in 
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FIG. 1: (Color online) Beam splitter array with each BS de- 
scribed by the unitary matrix U, Eq. |TJ. Particles start in modes 
with label (0), propagate from top to bottom through the ad- 
jacent rows of BS, and are eventually detected in modes with 
label (n), on exit. 

a coherent superposition of the two output modes, as 
described by the unitary 2x2 matrix 



U = 



1 



1 i 
i 1 



(1) 



2 



While the particles propagate through the array, we have 
to distinguish two inequivalent steps, as depicted in 
Fig. [TJ The direct sum Vi of the single BS matrices U 
describes odd time steps, where modes 1 and 2 enter 
one beam splitter, modes 3 and 4 enter the next and so 
forth. In the following step, modes 2 and 3 couple into 
the same BS, and so do modes 4 and 5 etc. Even steps 
are thus described by V2 = S~ViS + , with matrix elements 
Sf. = 5[(i - j + l) mo d2L,0] of S ± . The 2L x 2L matrix 



W(n) = 



( Vi V2 ) " , for n even 
(V 1 V 2 ) S 2 1 Vi, for n odd 



(2) 



describes the evolution for n time steps. The creation 



operator of a particle a^ ^ in mode j at time step is 

mapped onto a superposition of the creation operators 
at time step n, with 2L the maximum number of modes: 



2L 



a; (0) - 2> ]q (n)aJ n >. 

q=l 



(3) 



Given an N-particle initial state, distributed over 2L in- 
put modes, 
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Eq. |3} implies the evolution into the output state 

2L „ ( 2L V> 



i^„>=n4j E w H (n)a q n) 

j=i V r r U=i ) 



10). 



(5) 



In Q [5jl, |0) and rj are the vacuum state and the occu- 
pation number of mode j (such that Yl=i r j = N), respec- 
tively. 

III. MANY-PARTICLE OBSERVABLES 

Note that the number of terms in d5} grows exponen- 
tially with N, so that simulations of the dynamics of 
many-particle QWs are practically impossible for larger 
systems. However, expectation values of many exper- 
imental observables can be derived without complete 
knowledge of the full final state. One class of such ob- 
servables which reveals many -body interference are par- 
ticle number correlation functions a f ' n ^ . . . a + ^ n 'a^ n ' . . . a? n ' 
of m output modes ij. For m=l, this is the single-mode 
particle expectation number, ni = a[ai, simply given by 
the sum of 2L single-particle densities of a single-particle 

QW, 



2L 

(m) = £ |W] 



ki| 2 rk . 



(6) 



where Wy is the transition amplitude from mode k 
to mode i, and we left out the dependence of Wy on 
n, for simplicity of notation. Eq. (|6j holds for both, 
fermions and bosons, as well as for distinguishable par- 
ticles. Within these different particle-families, various 
distinct initial states are mapped on final states with 
the same mean particle number by the beam splitter ar- 
ray. Many-particle interference does not manifest on this 
level, and no information about the many-particle coher- 
ence of the state is obtained, as anticipated above. This 
is different when we correlate detection events in two 
different modes, i.e. m=2 (the upper/lower ±-sign refers 
to bosons/fermions): 

(•MW^'^Vfin) 
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= Yj l WkiWl i ±Wk j W li | 2 r k ri + Y |WkiW kj | 2 r k (rk - 1). 

(7) 
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The two sums in this equation represent two physi- 
cally distinct processes: in the first, the particles origi- 
nate from different modes. The two-particle amplitude 
WkiWij describes the simultaneous transition of one par- 
ticle from mode k to i, and of another one from 1 to j. Due 
to the indistinguishability of the particles, this amplitude 
interferes with that where k and 1 are exchanged. The 
second sum accounts for particles that originate from the 
same mode. The corresponding many-particle amplitude 
W k jW k j is invariant upon exchange of the input modes, 
and hence exhibits no many-particle interference. 

Eqs. l[6j ffj exhaust all possible correlations between 
two particles injected into the BS array. Only higher or- 
der correlations between a larger number of particles can 
exhibit genuine many-particle interference effects. Since 
the number of distinct correlations increases like the fac- 
torial of the particle number N, a general expression for 
the m-mode particle number correlation function is de- 
sirable, and can be shown to read 
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where ?i = (n — Si,q k )- The function fB/F(cr fl ) is given 

k=l 

by unity for bosons (B), and by the sign of the permu- 
tation in its argument for fermions (F). The outer, inco- 
herent sum in Eq. (JSj> runs over all ordered combinations 
q = (qi, ..,q m ) of input modes the particles may have 
originated from. Every single summand represents a 
physically distinct event. Each such event, in turn, is 
given by the coherent sum over all S(q) permutations cj„ 
of q, due to the particles' indistinguishability, where 0„j 
represents the j-th element of o^. Clearly, ||8j| reproduces 
J6) and |7lfor m=l and m=2, respectively, with the two 
sums in^7| corresponding to qi < q 2 and qi = q 2 in {§}. 
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Expression dsl permits to infer a large variety of 
many-particle observables, including the particle num- 
ber counting statistics. Since the latter are restricted to 
occupation numbers zero or one for fermions, we con- 
centrate on bosons in the following. 



IV. COUNTING STATISTICS 

The counting statistics of mode i, i.e. the probability 
P^(k) to find k particles in mode i, is defined by the 
linear relation 

N 1 i 

Qm := (a^ar > = £ P^W^-^ , < m < N, 

k=0 ^ 

(9) 

where the are determined via Eq. (pi. This relation 
can be inverted to yield P^(k), if all moments are 
known. Analogously the counting statistics P^^k^kj) 
of two modes i and j, i.e. the joint probability to find 
exactly ki particles in mode i and kj particles in mode j, 
is given in terms of 

rfi$ - / a +W m ' +W m i (n) m ' (n) m i s 

Virgin, •- \ a i a j a i a j / 

- y p(M) (k . k .)^ V_ (10) 

with < mi, mj < N. Thus all two-mode moments Qm|' mj 
are required for the two-mode counting statistics. Multi- 
mode counting statistics for more than two modes are 
obtained in strictly analogous fashion. 



A. Single-mode statistics 

With these tools, we can now fully characterize the dis- 
tribution of particles among the modes, and work out the 
observable consequences of many-particle interference: 
We start out with N=8 particles, launched in eight ad- 
jacent modes at the center of an array of length 2L=50, 
and initially in a product state (|4j. Fig.|2]shows the prob- 
ability for mode 25 to be occupied byk particles after 
n = 6 steps. The probability decreases with k, both for 
bosons and distinguishable particles. Since, however, 
for bosons all many-particle amplitudes that contribute 
to the coherent sum in ||8| are identical, they interfere con- 
structively. Consequently, the probability to find more 
than two particles in one mode decreases much slower 
for bosons than for distinguishable particles. For the 
latter, P' 25 ^(k) exhibits Poissonian decrease like 1/k!. In 
contrast, in the limit of many particles, P' 25 '(k) can be 
shown to approach exponential decrease in the bosonic 
case, which is apparent in Fig. [2] for low and intermedi- 
ate values of k. For k=N=8, p( 25 '(k) is equal to QSj_ 8 in 
(|9j. There, all N! amplitudes in the coherent sum in ||8| 
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FIG. 2: (Color online) Single-mode counting statistics P (25) (k) 
for mode i=25 after n = 6 steps on a logarithmic scale. Bosonic 
bunching enhances P (25 '(k) for bosons by several orders of mag- 
nitude (blue circles) as compared to the result for distinguish- 
able particles (red squares). The inset shows the distribution 
for small k. 

are necessarily equal. Thus, bosonic bunching 112311241 in- 
duces a factor N! =8! as compared to distinguishable par- 
ticles. For k < N, P< 25 >(k) is a sum of several Q^, what 
prevents such simple derivation of the bosonic value of 
P< 25 >(k). 

Further note that, since the number of constructively 
interfering amplitudes grows monotonically with the 
number of detected particles, the distribution in Fig. [2] 
is smooth (in contrast to the strong fluctuations in Fig. 
below), and, in particular, it is qualitatively independent 
of the considered output mode. 



B. Two-mode statistics 

In contrast to the above single-mode counting statis- 
tics, distinct many-particle amplitudes contribute to the 
coherent sum in d8l when several output modes are corre- 
lated, and give rise to characteristic many-particle inter- 
ference effects beyond the smooth bunching observed in 
Fig. [2] For these interference effects to fully develop, 
we need to propagate the many-particle state over a 
larger number of steps (n=20), such that all single par- 
ticle wave functions spread over the entire lattice, and 
mode-to-mode correlations can build up. Our figure 
of merit is then the conditional probability P^'(ki,kj) to 
find ki/j particles in mode i/j, with a total of m = kj + k 
particles in these modes. Indeed, as evident in Fig. [3] 
many-particle interference of indistinguishable particles 
has a marked impact on P^ N=g (ki,kj) as compared to 
the case of distinguishable ones. While the latter always 
implies a single-peaked dependence of Pj^ N=g (ki,kj) on 
the particle imbalance Ak = kj - kj, many-particle inter- 
ference manifests in strong modulations of P^ N=g (ki, kj) 
for bosons, with pronounced maxima for even values of 
Ak, i.e. for odd particle numbers in either one of both 
modes (remember that N=8). Furthermore, the many- 
particle interference condition is strongly affected by the 
specific choice of the output modes i and j, as obvious 
from a comparison of Figs. [3^) and b): an asymmetric 
choice of i and j with respect to the lattice center induces 
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FIG. 3: (Color online) Joint counting statistics for two different 
pairs of modes: Particle distribution P2°'(Ak = k ; - kj) among 
two modes after n=20 steps of N=8 particles, for lattice size 
2L=50, conditioned on finding m=4 or 8 particles, a) Joint 
probability distribution of the symmetric (with respect to the 
lattice center) modes i=19, j=32 to find m=4 particles (dashed; 
distinguishable: red squares, bosons: blue circles) and m=8 
particles (solid; distinguishable: red squares, bosons: blue cir- 
cles), b) Same as a), for asymmetric modes i=18, j=32. c) 
P™ N (Ak = kj - kj) of modes i=19, j=32 for N=4 (orange dia- 
monds), N=6 (green squares) and N=8 (black circles). 



single-peaked and essentially n-independent for distin- 
guishable particles, some of the interference structure 
of Fig. [3] prevails under the average, and, furthermore, 
shows considerable variation with n. This latter depen- 
dence is, similarly to the case considered in Fig. [3j re- 
duced for m < N. 
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FIG. 4: (Color online) Two-mode statistics as in Fig.[3j averaged 
over all pairs of modes, after a) n=19 steps and b) n=20 steps. 
Color code as in Fig. [3] Note that while the distribution is 
essentially invariant under changes of n for distinguishable 
particles, the correlations fluctuate with n, in the bosonic case. 



a strong asymmetry in F m (ki, kj). What's more, the in- 
terference signal sensitively depends on the total num- 
ber of particles N in the system, in a non-monotonic way, 
as apparent from Fig. [3] c), where the bosonic signal is 
shown for different N. This allows for a clear distinction 
of multi-particle interference contributions of different 
order. 

When not all particles contribute to the signal in modes 
i and j, i.e. for m = ki + kj < N, many of the events which 

contribute to Pj^ N (ki,kj) can be distinguished by the 
final positions of the remaining (N - m) particles and 
thus cannot interfere anymore, what leads to a smoother 
dependence of P^ N (kj,kj) on Ak. This is clearly ap- 
parent in Figs. |3]a),b), as well as from |[8}: for expecta- 
tion values that involve the maximum number N = 8 of 
creation/annihilation operators, the incoherent sum col- 
lapses to a single term while the coherent sum contains 
up to N! terms. For m < N, many more incoherently 
added terms contribute, and, thus, interference effects 
are averaged out. 

Let us finally stress that the observed many-particle 
interference effects even persist after an average of 

P^= N=8 (ki,kj) over all pairs (i,j): Fig. 4 shows that the 
bosonic signal is markedly different from the classical 
result, after a given number of steps as well as in its 
dependence on n. While the two-mode correlation is 



V. CONCLUSION 

In summary, many-particle interferences between N 
indistinguishable particles injected into the 2L input 
modes of a beam splitter array give rise to strong bunch- 
ing effects in the single mode counting statistics, as well 
as to a robust, pronounced and non-monotonic depen- 
dence of the two-mode correlation function on total 
particle number and particle imbalance between both 
modes. The many-particle interference correction in- 
duces large deviations from the classical (distinguishable 
particle) result. This allows to characterise the witnessed 
interference phenomenon in terms of the number of con- 
tributing particles - an information unaccessible, e.g., in 
few-particle quantum walks. This is of particular inter- 
est, e.g., in the many-particle transport across disordered 
media |25|, which is also covered by our central result, 
Eq. |8|, that actually allows to explore the full abundance 
of statistical observables of any non-interacting system 
with an arbitrary number of particles. 
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